

  filenamo = "radon_obs_mm.nc" 

  system("rm "+filenamo) 

  filo = addfile(filenamo,"c") 


  sitelist = asciiread("radon_sitename.dat",-1,"string")   
  siteeur  = asciiread("radon_siteeur.dat",-1,"string")   
  siteasia = asciiread("radon_sitechn.dat",-1,"string")   
  sitenam  = asciiread("radon_sitenam.dat",-1,"string")   
  siteplot = asciiread("radon_siteplot.dat",-1,"string")   

  ns = dimsizes(sitelist) 


  str = new(10,"string") 
  str = "" 

  rmm = new((/12,ns/),"float") 

  rmm!0 = "time" 
  rmm&time = ispan(1,12,1) 
  rmm!1 = "site" 
  rmm&site = ispan(1,ns,1) 
  rmm@long_name = "monthly mean radon conc." 
  rmm@units = "mBq/m3 STP" 
  rmm = 0.  

  keur = new((/ns/),"integer") 

  keur!0 = "site" 
  keur&site = ispan(1,ns,1) 
  keur@long_name  = "Site in Europe (Y:0,N:1)" 
  keur(:) = (/ siteeur(:) /) 

  kasia = keur 
  kasia@long_name = "Site in Asia (Y:0,N:1)" 
  kasia(:) = (/ siteasia(:) /) 

  knamr = keur 
  knamr@long_name = "Site in North America (Y:0,N:1)" 
  knamr(:) = (/ sitenamr(:) /) 

  kplot = keur 
  knamr@long_name = "Site for plot (Y:0,N:1)" 
  kplot(:) = (/ siteplot(:) /) 
 
  do is = 0,ns-1
     print(""+sitelist(is))  
     rmn = asciiread("./data/radon.obs."+sitelist(is)+".9999.dat",-1,"float") 
     rmm(:,is) = (/ rmn(:) /)  
  end do 

  filo->radon=rmm
  filo->site_EUR=keur
  filo->site_ASIA=kasia
  filo->site_NAMR=knamr
  filo->site_plot=kplot




